AD  747031 


AFCRL-72-0077 


( 

I 

Distribution  of  Water  Vapor  in  the  Stratosphere  as 
Determined  from  Balloon  Measurements  of  Atmospheric 
Emission  Spectra  in  the  24  -  29/Jim  Region 

by 

Aharon  Goldman,  David  G.  Mure  ray,  Frank  H.  Mure  ray 
Walter  J.  Williams,  James  N.  Brooks 


Department  of  Physics 
University  of  Denver 
Denver,  Colorado  80210 


Contract  No.  F19628-71-C-0171 
Project  No.  8692 


Scientific  Report  No.  1 
February,  1972 


D  D  Ov 


InUS? 


■\mnn 


Pi  AUG  22  1972 


Contract  Monitor:  Robert  A.  McClatchey 
Optical  Physics  Laboratory 


C 


Approved  for  public  release;  distribution  unlimited. 


Sponsored  by 

Advanced  Research  Projects  Agency 
ARPA  Order  No.  1366 


Monitored  by 

AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS  01730 


Reproduced  by 

NATIONAL  TECHNICAL 
INFORMATION  SERVICE 

U  S  Deportment  of  Commerce 
Springfield  VA  22151 


I 


£Z  w.  i 


.’■v' A  ■»  •$* 

I 

eP“ 


I?  n  l  i 

.  -  • 


rfll 


Program  code  no . . . .  IE50 

"C '  -  ” 

Effective  date  of  contract . .  26  March  1971 

Contract  expiration  date . . . . .  29  February  1972 

Principal  investigator  &  phone  no . Dr.  David  G.  Murcray/303  753-2627 


Project  Scientist  &  phone  no. . . .  JDr.  Robert  A.  McClatchey/617  861-3224 


1  MCESICm  lir  1 

OTIS 

Wfclti  Jtctki 

DOC 

Bnff  Sectiu 

D 

WMNMV.UEl 

q 

JIBUflCATlOH  ... 

V . 

BISTRISUTION/SVSIUSILIIT  COSES 

.  Olsl  "  AVAIL.  i;id/of  SfECMl 

Qualified  requestors  may  obtain  additional  copies  from 
the  Defense  Documentation  Center.  All  others  should 
apply  to  the  National  Technical  Information  Service. 


Unclassified 


Security  Classification  - - 


DOCUMENT  CONTROL  DATA  ■  R  &  D 

(Sacurity  claaallleatlon  of  tltla,  body  ol  abalract  and  tndexlnj  annotation  mu»l  ha  entered  whan  th^ovarall_ra£ortja_claaalllad2_ 


I.  originating  ACTIVITY  ({.jrporata  author) 

University  of  Denver 
Department  of  Physics 

Denver,  Colorado  80210  _ 


3.  REPORT  TITLE  _  _  .  _ _ _ 

DISTRIBUTION  OF  W  VTER  VAPOR  IN  THE  STRATOSPHERE  AS 
DETERMINED  FROM  BALLOON  MEASUREMENTS  OF  ATMOSPHERIC 
EMISSION  SPECTRA  IN  THE  24gm  -  29um  REGION _ . 


4  deSvSipti  ve  notes  (Typa  ol  rmport  and  Inclualve  dataa) 

Scientific  Interim  _ _ 


5-  AUTHOR(S)  (flrat  name,  mlddtm  Initial,  laat  name) 

Aharon  Goldman  Walter  J.  Williams 

David  G.  Murcray  James  N.  Brooks 

_ Frank  H.  Murcray _ _ _  ■ - 

t  REPORT  date: -  7a.  TOTAL  NO.  OF  PAGES  |?b.  NO.  OF  REFS 

February  1972  _  38 


6a.  CONTRACT  OR  GRANT  NO. 

F19628-71-C-0171  ARPA  Order  No.  1366 

b,  PROJEC  T  NO. 

Protect,  Task,  Work  Unit  Nos. 
c  8692  n n/a 

DoD  Element  6230 ID 


9a.  ORIGINATOR'S  REPORT  NUMBER(S) 


Scientific  Report  No.  1 


96.  other  REPORT  NOIS)  (Any  other  numb*1  a  that  may  be  aaalgnad 
thla  report  \ 


AFCRL-72-0077 


10.  DISTRIBUTION  STATEMENT 


A  -  Approved  for  public  release;  distribution  unlimited. 


II.  SUPPLEMENTARY  NOTES 


12.  SPONSORING  •••'.!  TAR"  ACTIVITY 


13  ABSTRACT 


This  research  was  supported  by  Ail.  Force  Cambridge  Research 
the  Advance  Research  Projects  Laboratories  (OP) 

AgencY- _ ' _ L  U  G.  Han  scorn  Field  _ 

Bedford,  Massachusetts  01730 

The  stratospheric  water  vapor  mixing  ratio  altitude  profile  has 

been  derived  from  spectral  observations  of  the  downward  night  emission 

from  the  pure  rotation  water  vapor  lines  in  the  24  -  29^m  region  of  the 

spectrum.  The  data  were  obtained  during  two  balloon  flights,  made  on 

22  February  1971  and  on  29  June  1971,  using  a  balloon-borne  spectral 

radiometer  with  ~2  cm’1  resolution.  The  observed  radiances  have  been 

fitted  to  line  by  line  -  layer  by  layer  radiance  calculations,  from  which  the 

water  vapor  mixing  ratio  between  10  and  30  km  has  been  derived.  The 

resulting  mixing  ratio  altitude  profiles  from  both  flights  show  a  broad 

minimum  around  15  km  and  a  broad  maximum  around  25  km, with  a  range 

7  6 

of  values  of  6  x  10  g/g  to  4  x  10  g/g. 


.1473 


Unclassified _ 

Security  Classification 


Security  CUeelfication 


Infrared  Emission  Spectra 
Line  by  Line  Model 
H_0  Mixing  Ratio 

£4 


AFC  RL- 72 -0077 


Distribution  of  Water  Vapor  in  the  Stratosphere  as 
Determined  from  Balloon  Measurements  of  Atmospheric 
Emission  Spectra  in  the  24  -  29fim  Region 

by 

Aharon  Goldman*,  David  G.  Murcray,  Frank  H.  Murcray 
Walter  J.  Williams,  James  N.  Brooks 


Department  of  Physics 
University  of  Denver 
Denver,  Colorado  80210 

Contract  No.  F19628-71-C-0171 
Project  No.  8692 


Scientific  Report  No.  1 
February,  1972 

*On  leave  from  the  Department  of  Environmental  Sciences,  Tel- Aviv 

University,  Tel-Aviv,  Israel 


Contract  Monitor:  Robert  A.  McClatchey 
Optical  Physics  Laboratory 


Approved  for  public  release;  distribution  unlimited. 


Sponsored  by 

Advanced  Research  Projects  Agency 
ARPA  Order  No.  1366 
Monitored  by 

AIR  FORCE  CAMBRIDGE  RESEARCH  LABORATORIES 
AIR  FORCE  SYSTEMS  COMMAND 
UNITED  STATES  AIR  FORCE 
BEDFORD,  MASSACHUSETTS  01730 


1 


ABSTRACT 

The  stratospheric  water  vapor  mixing  ratio  altitude  profile  has 

been  derived  from  spectral  observations  of  the  downward  night  emission 

from  the  pure  rotation  water  vapor  lines  in  the  24  -  2%m  region  of  the 

spectrum.  The  data  were  obtained  during  two  balloon  flights,  made  on 

22  February  1971  and  on  29  June  1971,  using  a  balloon-borne  spectral 

radiometer  with  ~2  cm  *  resolution.  The  observed  radiances  have  been 

fitted  to  line  by  line  -  layer  by  layer  radiance  calculations,  from  which  the 

water  vapor  mixing  ratio  between  10  and  30  km  has  been  derived.  The 

resulting  mixing  ratio  altitude  profiles  from  both  flights  show  a  broad 

minimum  around  15  km  and  a  broad  maximum  around  25  km  with  a  range 

7  6 

of  values  of  6  x  10  g/g  to  4  x  10  g/g. 
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I.  INTRODUCTION 

A  number  oi  spectrrjscopic  derivations  of  the  water  vapor  vertical 
distribution  in  the  atmosphere  have  been  made.  The  earlier  spectroscopic 
measurements  were  made  from  observations  of  the  absorption  of  the  solar 
infrared  spectrum  from  aircraft^  ^  and  balloons.  More  recently, 

water  vapor  amounts  have  been  inferred  from  infrared  radiometric  mea¬ 
surements  of  the  water  vapor  emission  spectrum.  These  include  data 

from  balloon-borne  radiometers  used  to  observe  the  infrared  downward 

(20-24) 

and  upward  emission  and  satellite -borne  radiometers  used  to  ob- 

..  ,  .  .  (25-31)  .  ..  ,  .  .  (32,33) 

serve  the  upward  emission  and  the  limb  emission. 

H^O  mixing  ratios  derived  from  these  data  show  considerable  di¬ 
vergence,  especially  for  the  middle  stratosphere.  The  results  for  the 
lower  stratosphere  are  less  contradictory  and  most  of  the  spectroscopic 

derivations  yield  a  "dry"  lower  stratosphere,  with  mixing  ratios  close  to 

-  6 

2  x  10  g/g.  Almost  all  of  the  frost  point  measurements  support  these 
(34-37) 

results.  However,  despite  the  numerous  experiments,  the  distri¬ 

bution  of  water  vapor  in  the  middle  stratosphere  is  rather  uncertain. 

Some  results  indicate  that  the  mixing  ratio  remains  close  to  the  "dry"  value 
of  the  lower  stratosphere,  while  others  indicate  that  the  mixing  ratio  in¬ 
creases  again  at  the  middle  stratosphere.  Comparison  of  the  different 
results  leads  to  uncertainties  of  one  to  two  orders  of  magnitude. 

It  is  evident  that  more  accurate  determinations  of  H  O  concentrations 

Cd 

in  the  middle  stratosphere  are  needed.  This  is  difficult  to  obtain  with  the 

absorption  spectra  method.  The  basic  limitation  of  this  method  is  that 

such  dat^  even  with  high  resolution  of  ~0.  25  cm  *  as  obtained  during  re- 

08  39) 

cent  balloon  flights,  '  '  cannot  be  used  to  meaningfully  measure  spec¬ 

trally  degraded  emissivities  which  are  smaller  than  a  few  per  cent.  On 

the  other  hand,  a  sensitive  cooled  spectral  radiometer  can  measure  spec- 

_3 

trally  degraded  emissivities  of  <10  .  Thus,  if  the  spectral  radiometer 
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can  resolve  line  groups,  the  spectral  emission  method  offers  significant 
advantages  for  the  determination  of  small  gas  amounts  in  the  upper  at¬ 
mosphere. 

In  the  present  study,  the  infrared  spectral  emission  method  is  ap¬ 
plied  to  the  derivation  of  the  H20  distribution  in  the  stratosphere  using 
the  spectral  radiance  data  obtained  during  two  recent  balloon  flights,  made 
on  February  22,  1971,  and  June  29,  1971,  from  Holloman  Air  Force  Base, 
New  Mexico.  During  these  flights,  spectral  observations  of  atmospheric 
downward  night  emission  were  obtained  in  the  16  -  30p,m  region  of  the 
spectrum.  Spectrsl  scans  were  made  from  a  number  of  altitudes  up  to 
30  km,  with  a  balloon-borne  grating  radiometer  system.  The  grating 
spectrometer  was  cooled  with  liquid  nitrogen  to  reduce  the  thermal  back¬ 
ground,  thus  increasing  the  sensitivity  of  the  liquid  helium  cooled  Cu:Ge 
detector.  The  resolution  obtained  is  better  than  2  cm"1  near  25ym. 

The  16  -  30^m  region  of  the  atmospheric  emission  spectrum  con¬ 
tains  line  features  due  to  CO?,  NzO,  HNC>3  and  H20,  as  has  been  verified 
from  the  high  resolution  solar  absorption  spectra  observed  in  this  region 

(39) 

wiih  a  balloon-borne  spectrometer.  In  the  long  wavelength  interval  of 

this  region,  between  24^tm,  and  30y,m,  the  effect  of  CO  ,  N  O  and  HNO  is 

2  2  3 

negligible,  and  the  observed  emission  is  mainly  due  to  pure  rotation  H  O 

£ 

linf  s. 

The  individual  H20  line  parameters  (positions,  intensities,  half¬ 
widths  and  energy  levels)  are  known  for  this  spectral  region,  hence  it  is 
possible  to  derive  the  atmospheric  vertical  distribution  of  H  O  on  a  line 
by  line  basis.  The  measurements  were  obtained  from  many  altitudes,  so 
that  the  atmosphere  can  be  divided  into  a  number  of  layers  and  the  line  by 
line  technique  can  be  applied  to  a  layer  by  layer  inferrenco  of  the  H  O 
profile.  Such  a  method  of  analysis  is  described  in  Section  II.  The  resul¬ 
ting  profiles  from  the  two  balloon  flights  are  described  in  Section  III,  fol¬ 
lowed  by  a  discussion  of  the  technique  in  Section  IV  and  the  conclusions  in 
Section  V. 
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II.  METHOD  OF  ANALYSIS 


The  radiometric  measurements  are  made  from  different  altitudes 

h  h . h  at  a  constant  zenith  angle  of  45  .  The  atmosphere  is  as- 

12  N 

sumed  to  be  stratified  into  layers  defined  by  actual  altitudes  of  observa¬ 
tion.  The  analysis  starts  at  the  highest  altitude  of  observation  h  ,  at  the 
smallest  radiance  values,  and  then  proceeds  layer  by  layer,  to  the  lower 
altitudes  with  increasing  radiance  values.  The  residual  atmosphere  above 

h  is  considered  as  one  homogenous  layer,  so  that  the  infinite  resolution 
1  - 1 

radiance  at  wavenumber  \(cm  )  to  be  observed  at  altitude  h  ,  R(v, h^),  is 

A  (40) 

a  one  layer  emission  according  to 


R(v,h1)  =  e(v,  1)  B(m  1),  (1) 

where  B(v,  1)  is  the  black  body  radiance  at  temperature  0 of  layer  1,  and 
e(v,  1)  is  the  infinite  resolution  emissivity  of  this  layer.  The  atmosphere 
below  hj  is  divided  into  N-l  homogenous  layers,  defined  by  h^,  h^,  .  .  .  ,  h^ 
The  emission  from  the  first  layer  is  transmitted  through  the  second  layer 
according  to  t(v,2)  R(v,  h^,  where  t(v,  2)  is  the  infinite  resolution  trans¬ 
mittance  of  the  second  layer.  In  addition,  the  second  layer  emits  radi¬ 
ance  according  to  e(n  2)  B(v  ,  2).  Thus,  the  infinite  resolution  radiance 
observed  at  h2,  R(v,h2),  is  given  by 

R(v,h  )  =  t(v,  2)  R(v,h1)  +  e(v,  2)  B(v,  2).  (2) 

This  process  proceeds  through  all  the  layers  below  h^  according  to 

R(v,  h  )  =  t(v,  n)  R(v,  h  )  +  e(v,n)  B(v,  n),  (3) 

n  n-i 

for  n  =  2, .  . .  ,  N. 

Eqs.  (1)  -  (3)  represent  successive  equations  for  R(\*h^),  where  each  step 

below  hj  involves  two  consecutive  layers.  (See  Appendix) 

Assuming  no  scattering,  the  relation  t(v)  =  l-e(v)  will  hold  for 

thermodynamical  equilibrium  and  this  can  be  substituted  into  Eqs.  (1)  -  (3). 

T(uh  )  is  a  function  of  the  amount  of  water  vapor  U(n)  in  the  layer  n  between 
n 
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altitudes  and  .  It  seems  now  that  once  a  theoretical  model  is  de¬ 
veloped  for  t(v, h^),  Eqs.  (1)  -  (3)  can  be  solved  successively  for  U(n), 
n  =  1, .  . . ,  N,  thus  yielding  the  desired  H^O  profile.  However,  the  obser¬ 
ved  spectral  radiance,  R(y,hn),  is  not  at  infinite  resolution,  but  is 
rather  a  convolution  of  the  infinite  resolution  radiance  R(p,hn),  with 
the  spectrometer  slit  function  g(V-V')  according  to^°^ 

R(v’V  = /yta  R(v'V  g(v-v')dv'/J^  g(v-v')dv',  (4) 

where  2a(cm_1)  is  the  spectral  width  of  the  slit  function. 

The  integrated  radiance  though,  is  independent  of  the  slit  func¬ 
tion^0)  so  that: 

Jvj  R(v«hn>dv  =  R(v,hn)dH  (5) 

where  v  ^  and  are  the  limits  of  the  spectral  region  under  consideration. 
Thus,  the  integrated  radiance  is  a  suitable  parameter  for  a  fitting  of  theo¬ 
retical  integrated  radiance,  Jtyv.h  )dv,  to  the  experimental  values, 
J>R(v.hn)dv,  by  successive  solutions  of  the  equations: 

V  v 

Ftu(ft)]  =  JVi  R(v,hn)dv  -  ^  R(v,hn)dv=  0,  (6) 

for  U(n),  n  =  1, .  . .  ,  N,  where  R{v,  h^)  is  given  by  Eqs.  (1)  -  (3). 

It  is  not  necessary  here  to  use  the  entire  band,  but  only  to  pick  out 
an  interval  [v^,  with  a  representative  isolated  group  of  lines.  It  is  im¬ 
portant  that  the  spectral  region  chosen  will  be  isolated,  so  that  the  experi¬ 
mental  integrated  radiance  can  be  accurately  determined. 

The  above  determination  of  U(n)  requires  a  suitable  model  for  the 
calculation  of  T(v,n).  The  most  accurate  transmittance  is  obtained  by  a 
line  by  line  calculation.  The  monochromatic  transmittance  t(v,  n)  is  given 
by(4l) 


T(v,n)  =  exp[-k(v,  n)  U(n)], 


(7) 
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where  k(v,n)  is  the  absorption  coefficient,  which  is  a  sum  of  the  absorption 

coefficients  of  all  spectral  lines, 

k(v,  n)  =  ^  k.(v,  n),  (8) 

i 

which  might  contribute  to  the  transmittance  at  v. 

The  atmospheric  infrared  line  shapes  are  given  by  the  Voigt  pro¬ 
file,  (4°'41)  which  is  a  convolution  of  the  Lorentz  and  Doppler  shapes.  The 

Doppler  halfwidth,  a^,  is  given  by 

a  (cm-1)  =  3.  58  x  10"7  v  (cm-1)  [6(K)/M]1/2,  (9) 

D  u 

where  vQ  is  the  line  position,  T  the  temperature  and  M  the  molecular 

weight.  At  6  =  225K,  v=  350  cm-1,  the  Doppler  half  width  is  _ 

q  =  0.  0004  cm-1.  Typical  collision  halfwidths  are  a  ~0.  08  cm  atm 
D 

Thus,  a  ~a  only  at  about  5  mb  (~36  km)  and  Doppler  effects  can  be  ne- 
'  D  L 

glected  at  the  altitudes  of  interest  (10-30  km).  Thus, 


k(  v,  n) 


S°.(n) 


a.(n) 

_ 1  9 

2  2 

(v-v  .)  +  a  .(n) 

oi  i 


(10) 


where  v  ,  a. ,  and  S°.  are  the  individual  line  centers,  line  halfwidths  and 
oi  i  i 

line  intensities  respectively.  The  individual  line  parameters  for  the 

region  of  interest  have  been  compiled  under  the  AFCRL  Atmospheric  Trans¬ 
mittance  Program.  (42)  In  the  atmosphere,  a.  and  S°.  will  vary  from  one 
layer  to  another  as  a  result  of  temperature  and  pressure  variations.  The 
line  intensities  are  assumed  to  depend  on  the  temperature  ©according 


S°(0 )  =  S°(6o)  [0o/e]3/2  exp[-1.439E"(l/0 


l/eo)][l-exP(-l.  439vq/0)] 


/[l-exp(-l.  439vo/0q)], 


(ID 
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where  E"  (cm  )  is  the  lower  state  energy,  0  ,  0  are  in  Kelvins,  and 

-1  ° 

Vo  ^cm  )  1®  *-he  Hne  position.  The  induced  emission  term,  normally  ne¬ 
glected  in  the  near  infrared  region,  is  included  here,  since  it  becomes 
significant  at  these  longer  wavelengths.  At  0  =  225K  and  0  =  300K, 

l-exp(-l.  439  v  70)  ~  0.96  and  ~  0.91  respectively  even  for  v  =  500  cm"\ 
o  o 

The  pressure  dependence  of  a  in  an  homogenous  layer  is  a  =  a°P 
where  P  is  the  effective  pressure  for  H  O  -  N  mixture,  given  by  ^3"45) 

Pe  =  Pt+4p'  (12) 

where  P  is  the  total  pressure  of  the  gas  mixture  and  p  is  the  partial 
pressure  of  the  absorbing  gas.  For  atmospheric  ^O,  p  is  small  so  that 
Pg  ~P^.  The  temperature  dependence  of  a  is  assumed  to  be^4^ 

<*°(0)  =  ot°(  0  )  (0  /0)°-62.  (13) 

o  o 


Atmospheric  temperature  and  pressure  profiles  can  be  obtained 

from  sondes  data  taken  during  the  radiance  measurements.  The  Curtis- 

(41)  _ 

Godson  approximation  is  then  used  to  derive  average  temperatures  0 

and  pressure  P  for  each  homogenous  layer.  A  constant  H  O  mixing  ratio 

11  (41)  2 

is  assumed  within  each  layer,  so  that 


P  =  0.  5(P  +P  ), 
n  n  n- 1 


(14) 


where  P  is  the  atmospheric  pressure  measured  at  altitude  h  .  A  similar 
n  _  n 

relation  is  assumed  for  the  temperature  0^.  These  expressions  for 
average  pressures  and  temperatures  are  used  for  n  >  2  while  for  n  =  1 
Pj  is  taken  as  P^/2  and  0^  as  0^. 

With  the  assumptions  made  above  for  the  transmittance  calculation 

one  can  proceed  to  the  solution  of  Eq.  (6).  It  turns  out  that  the  Newton- 

(47) 

Raphson  iterative  method  is  very  efficient  for  the  solution  for  U(n). 

The  iterations  are  made  according  to 

•  i  1  • 

U1  (n)  =  U1  (n)  -  F[U(n)]  /  F'  [U(n)], 


(15) 
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where  the  prime  denotes  the  derivative  with  respect  to  U(n),  and  U  (n) 
is  an  approximation  for  the  zero  of  F[U(n)].  The  derivative  is  given  by 

v- 

F‘[U(n)]  =  -kK  n)  exp[-k(v,n)  U(n)]  R(v,hn_1)dy 
v 

+  f  k(v,  n)  exp[-k(v,n)  U(n)]  B(v,  n)dv,  (lo) 

V1 

for  n  >  1,  while  for  n  =  1  the  first  term  is  omitted.  The  absorption  co¬ 
efficients  k(v,  n)  have  to  be  computed  only  once  for  each  layer,  and  are 

not  changed  during  the  iterations. 

A  computer  program  has  been  developed  to  derive  the  HzO  profile 
according  to  the  method  outlined  above.  The  computati  on  of  k(v,  n)  has  been 

made  by  subroutines  previously  developed  by  T.  G.  Kyle.  The 

program  can  also  produce  a  degraded  emission  spectra  for  a  given 
profile.  Thus,  once  values  for  U(n)  are  derived,  the  calculated  degraded 
spectra  can  be  compared  with  the  measured  degraded  spectra,  for  any 

altitude  h  . 

n 

Despite  the  apparent  simplicity  of  the  equations  involved,  large 

high  speed  computers  are  needed  to  make  such  computations  practical.  A 

10  layer  calculation  for  a  15  cm  1  interval  with  4  iterations  per  layer  takes 

~60  seconds  on  the  CDC  7600  at  NCAR,  Boulder.  This  is  mainly  due  to 

the  fine  net  needed  for  the  calculation  of  the  absorption  coefficients  and  due 

to  the  fact  that  each  step  in  the  layer  by  layer  calculation  involves  the 

radiances  from  two  layers  (except  for  the  first  layer).  A  fine  net  interval 

of  0.  001  cm"1  was  needed  for  the  upper  layers  due  to  the  small  atmospheric 

«  (49) 
pressures. 

F[U(n)]  is  a  smooth  function  of  U(n),  so  that  with  a  reasonable 
guess  for  an  initial  approximation,  the  single  zero  of  F  can  be  found  to 
4  -  5  significant  figures  withing  3  to  5  iterations.  However,  an  improper^ 
value  for  the  initial  approximation  uV)  may  yieU  a  negative  value  for  U 
(n)  so  that  a  new  uV)  has  to  be  tried.  The  same  iteration  procedure  is 


repeated  for  each  new  layer,  until  the  complete  H^O  profile  is  derived. 
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HI.  RESULTS 


Selected  records  of  the  emission  spectra  obtained  in  the  16  to 

30pm  region  during  the  balloon  flight  made  22  February  1971  are  shown  in 

-7  -4  -2 

Fig.  1.  The  radiance  scale  for  each  record  is  from  10  to  10  w  cm 
sr"*p,"*  and  each  record  is  offset  from  the  next  by  two  decades.  The  zenith 
angle  was  45°  for  all  records.  The  times,  altitudes,  pressures  and  tem¬ 
peratures  are  shown  in  Table  1.  The  results  of  the  29  June  1971  flight 
show  similar  spectral  patterns,  but  with  somewhat  different  values  of 
radiances  as  a  function  of  altitude.  Temperature  and  pressure  data  for  the 
two  flights  were  obtained  from  rawinsonde  ascents  made  at  the  same  time 
and  place  as  the  primary  flights.  The  rawinsonde  data  obtained  during  the 
February  flight  are  not  complete,  so  that  it  was  necessary  to  supplement 
the  data  on  the  basis  of  rawinsonde  data  obtained  at  Holloman  on  different 
dates,  and  on  the  basis  of  the  U.  S.  Standard  Atmosphere.  Complete 

rawinsonde  pressure  and  temperature  data  were  obtained  during  the  June 
flight  and  these  values  were  used  for  the  data  analysis. 

Examination  of  Fig.  1  shows  that  for  wavelengths  larger  than  24pm 
the  radiance  is  dominated  by  I^O  pure  rotation  lines,  degraded  to  several 
line  groups  by  the  spectrometer  slit  function.  The  line  group  near  25pm 
is  particularly  well  isolated,  and  can  be  used  to  fit  JR(\,  h^JdX.  to  jR(\,hn)d\ 
as  described  in  Section  II.  The  spectral  region  of  this  group  extends  from 
389  to  405  cm"1,  and  contains  9  strong  H20  lines  and  18  much  weaker 
lines.  ^^Each  one  of  the  strong  lines  is  overlapped  to  some  extent  by  one 
or  more  of  the  other  strong  lines,  and  the  measured  radiance  shows  how 
these  are  degraded  to  3  distinct  peaks.  The  measured  integrated  radiance 
values  of  the  25pm  line  group  are  given  in  the  last  column  of  Table  I. 

Starting  with  one  layer  above  float  altitude,  U(l)  has  been  derived 
from  the  25pm  line  group,  yielding  a  mixing  ratio  of  5.  8  x  10  g/g  for 
the  February  flight.  This  value  of  U(l)  has  been  used  to  compute  a 
degraded  spectrum  for  the  20  -  30pm  region,  shown  in  Fig.  2.  It  is  seen 
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that  while  near  20(j,m  the  resolution  of  the  calculated  spectrum  is  higher  than 
that  of  the  experimental  spectrum,  the  opposite  is  the  case  near  30ym. 

This  is  due  to  the  fact  that  a  constant  wavenumber  slit  function  (in  this 
case  a  Gaussian  function  with  a  halfwidth  of  0.  90  cm’*)  has  been  used  for 
the  calculated  spectrum,  while  the  experimental  grating  spectra  has 
roughly  a  constant  wavelength  resolution.  The  integrated  radiances  are  not 
effected  by  the  slit  function  and  it  is  seen  that  even  though  U(l)  was  derived 
from  the  25(j,m  line  group  only,  the  agreement  between  the  calculated  and 
experimental  radiance  is  very  good  through  out  the  whole  region. 

Similar  comparisons  for  several  layer  calculations  over  short  inter¬ 
vals  within  the  24  -  29(jm  region  show  the  same  agreement.  In  particular, 

10  layer  calculations  for  both  the  25|im  line  group  and  the  26/im  line  group, 
yield  very  close  profiles.  Thus,  it  was  concluded  that  only  the  25(j,m  line 
group  will  be  used  for  the  derivation  of  the  water  vapor  profile.  This 
reduces  considerably  the  computer  time  and  allows  the  use  of  15  atmos¬ 
pheric  layers  in  the  computation  without  exceeding  1  minute  of  computing 
time  of  the  CDC  7600. 

The  results  for  U(n)  from  a  15  layer  calculation,  where  the  layers 
are  defined  by  the  records  in  Fig.  1  for  the  February  flight  are  shown  in 
Fig.  3.  The  histogram  representation  shows  the  thickness  of  the  layers 
and  stresses  the  fact  that  a  constant  mixing  ratio  is  assumed  within  each 
layer.  The  calculated  and  experimental  radiance  at  the  lowest  altitude, 
h15  =  8*  ^  km’  are  shown  in  Fi8-  4*  Results  from  a  15  layer  calculation 
for  the  June  flight  are  also  shown  in  Fig.  3. 

The  number  of  layers  was  limited  to  15  even  though  more  records 
were  obtained  at  intermediate  altitudes.  The  resulting  altitude  resolution 
shown  in  Fig.  3  is  about  2  km  near  float  altitude,  gradually  increasing  to 
~0.  5  km  near  10  km.  This  variable  altitude  resolution  has  been  chosen  on 
the  basis  of  the  diminishing  rate  of  the  radiance  decrease  with  altitude,  as 
evident  in  Fig.  1  and  Table  I.  The  altitude  resolution  will  increase  by  in¬ 
creasing  the  spectral  interval.  A  test  for  the  validity  of  the  February  pro- 


11 


file  shown  in  Fig.  3  was  made  by  deriving  a  10  layer  profile  in  which  some 
of  the  layers  from  the  15  layer  calculation  were  doubled  in  thickness.  The 
results  are  shown  by  a  broken  line  on  Fig.  3,  and  it  can  be  seen  that  these 
actually  represent  some  average  of  the  15  layer  profile. 

Some  data  are  missing  between  11  and  18  km  on  the  June  flight  so 
that  some  of  the  altitude  interval  sizes  are  not  optimum  whithin  this  alti¬ 
tude  range. 

Fig.  3  shows  that  for  February.  1971,  the  water  vapor  mixing  ratio 
above  29.  3  km  is  about  6  x  10-6  g/g  and  for  June,  1971,  about  4  x  10  g/g. 
The  present  method  does  not  infer  the  distribution  above  float  altitude  so 
that  these  values  represent  the  total  amount  of  water  vapor  above  float  if 
the  water  vapor  were  uniformly  mixed.  Both  the  February  and  June  pro¬ 
files  show  a  broad  minimum  extending  from  12  to  18  km,  with  a  mixing 
ratio  of  <  1  x  10  ^  g/g.  The  mixing  ratio  increases  above  this  altitude, 
forming  a  broad  maximum  near  24  or  25  km.  The  mixing  ratio  may  de¬ 
crease  again  above  25  km,  but  must  increase  at  some  higher  altitude  to 
account  for  the  observed  radiance  values  at  float  altitude.  The  minimum 
near  15  km  and  the  maximum  near  25  km  are  distinct  features  of  both 
flights,  even  though  the  minimum  and  the  maximum  are  lower  on  the 

June  flight. 
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IV.  DISCUSSION 


The  above  analysis  of  the  spectral  radiance  data  should  be  examined 

in  detail  to  determine  the  possible  sources  of  error  and  their  magnitude. 

Both  the  theoretical  parameters  and  the  experimental  data  may  contribute 

errors.  The  theoretical  parameters  are  discussed  in  detail  in  this  section 

and  wherever  possible,  the  experimental  parameters  are  also  considered. 

The  present  formulation  uses  an  isotropic  steady  state  model, 

which  does  not  take  into  account  the  time  and  position  variations  of  H  O 

2 

radiance  values  from  different  layers.  Once  a  U(n)  is  derived,  the  same 

U(n)  is  used  to  derive  U(n+1),  even  though  the  measurements  at  h  at  h 

n  n+1 

are  not  simultaneous  in  time  and  location;  they  are  made  during  ascent  or 
descent  of  the  balloon,  which  will  drift  a  number  of  miles  during  a  flight. 

The  present  model  applies  to  spectral  regions  which  are  dominated 
by  lines  centers  and  not  by  wings.  Emissions  from  clouds  or  particle  lay¬ 
ers  are  not  taken  into  account,  thus  assuming  a  "clear"  atmosphere.  This 
avoids  the  uncertainty  in  the  absorption  coefficient  of  the  continuum  and  the 
unsatisfactory  representation  of  the  line  wings  by  the  Lorentz  profile. 

Within  these  limitations,  the  iteration  procedure  can  fit  U(n)  to 
satisfy  Eq.  (6)  to  any  arbitrary  accuracy;  but  the  significant  digits  in  U(n) 

are  limited  by  the  accuracy  of  the  experimental  values  R(v,h  ),  of  the  in- 

n 

dividual  line  parameters  and  of  the  atmospheric  temperatures  and  pres¬ 
sures.  Line  position*  are  known  to  high  precision  (6  to  7  significant  fig¬ 
ures)  but  line  intensities  and  line  halfwidth  are  not  known  to  better  than 
10  to  20%.  This  obviously  limits  the  accuracy  of  the  transmittance  determ¬ 
ination.  The  uncertainties  in  atmospheric  temperature  and  pressure  data 
effect  both  the  transmittance  and  the  radiance  values.  Trial  runs  with  a 
number  of  atmospheric,  temperature  profiles  showed  that  errors  of  3  -  4K 

will  not  effect  the  derived  U(n)  by  more  than  5%.  The  error  in  U(n)  due  to 
c 

errors  in  S  and  ot  can  then  be  inferred  from  the  square  root  approximation 
for  a  single  line,  [S  aU(n)]  ^  .  Trial  runs  for  a  few  layers  with  dif- 
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ferent  values  for  S°and  a  have  verified  this  estimate  of  error.  Since  10% 
error  is  typical  for  S°and  a,  the  square  root  relation  may  increase  the  as¬ 
sociated  error  in  U(n)  to  40%. 

There  are  two  effects  these  errors  can  have  on  the  water  vapor 
mixing  ratio.  They  can  influence  the  absolute  accuracy  of  the  derived 
values  of  U(n)  and  of  the  relative  values  of  U(n)  at  different  altitudes.  The 
absolute  values  have  a  greater  uncertainty  than  the  relative  values.  This 
is  important  to  stress  because  the  uniqueness  of  the  present  altitude  mix¬ 
ing  ratio  profiles  is  in  their  relative  values  between  10  and  30  km.  The 
resultant  absolute  error, based  on  the  above  estimates  of  10%  for  S°  and  a 
5%  for  atmospheric  temperature  and  allowing  10  to  15%  error  in  the  exper¬ 
imental  measurement  of  R,is  up  to  60%  of  the  absolute  values  of  U(n).  The 
relative  values  seem  to  be  accurate  within  5  to  10%. 

In  evaluating  possible  errors,  one  additional  consideration  must  be 
made.  Previous  studies  have  shown  that  stratospheric  water  vapor  mea¬ 
surements  may  be  subject  to  errors  due  to  local  moisture  contamination. 
Most  of  these  studies  involve  measurements  of  the  water  vapor  content  in 
the  immediate  vicinity  of  the  instrument  package.  The  contamination  must 
occur  over  a  significant  path  length  before  it  is  important  in  an  infrared 
emission  study.  This  path  can  be  considered  in  two  parts,  that  which  is 
internal  to  the  instrument  and  that  which  is  external  to  it,  but  still  influ¬ 
enced  by  the  balloon  train  dynamics. 

Laboratory  studies  have  established  that  there  is  no  measurable 
water  vapor  contamination  within  the  spectral  radiometer.  This  is  due  to 
the  low  operating  temperature  of  the  radiometer  combined  with  the  nitrogen 
venting  system.  The  external  optical  path  is  relatively  short  compared 
with  the  path  from  the  atmospheric  layers  being  considered  and  would  re¬ 
quire  considerable  contamination  to  influence  the  results  significantly. 
There  is  good  evidence  to  indicate  that  this  amount  of  contamination  is  not 
present.  Water  vapor  emission  data  associated  with  the  flights  analyzed 
here  show  that  when  measurements  made  during  ascent  were  compared 
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with  those  made  during  descent,  the  values  were  identical  at  similar  alti¬ 
tudes.  Such  measurements  preclude  the  possibility  of  any  significant 
amount  of  water  being  swept  down  from  the  balloon  vehicle.  In  addition, 
similar  data  collected  over  an  hour  or  more  at  float  altitude  showed  no 
significant  increase  or  decrease  in  the  amount  of  vapor  in  the  emitting 
path. 
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V.  CONCLUSIONS 


On  the  basis  of  these  results,  some  of  the  fine  structure  in  the 
inferred  H^O  profiles  might  be  artificial,  but  the  gross  vertical  structure 
with  the  broad  maximum  near  25  km  seems  to  be  real.  The  maximum 
indicates  a  relatively  "wet"  atmospheric  layer  near  25  km.  The  origin  of 
such  a  layer  is  not  obvious  and  might  be  related  to  the  photochemical  pro¬ 
cesses  in  the  O  ,  HNO,,  NO_  layer.  Comparison  of  the  profiles  from  the 

J  J  M 

22  February  1971  and  29  June  1971  flights  indicate  that  a  seasonal  vari¬ 
ation  might  exist.  During  summer,  both  the  minimum  near  15  km  and  the 
maximum  near  25  km  are  smaller  than  that  during  the  winter.  Additional 
measurements  are  obviously  required  to  establish  seasonal  variations  in 
the  water  vapor  profile. 

The  procedure  developed  in  the  present  study  can  be  applied  to 
other  types  of  layer  by  layer  radiance  calculations.  The  inference  of 
H^O  profiles  from  upward  emission  measurements  from  balloons  or  satel¬ 
lites  can  be  accomplished  with  the  same  procedure  used  here,  only  rever¬ 
sed  in  the  order  of  altitudes,  and  the  earth  emission  term  has  to  be  inclu¬ 
ded.  Limb  radiance  scans  from  satellites  can  be  treated  in  a  similar  man¬ 
ner  even  though  the  geometry  of  the  experiment  is  different.  The  technique 
is  not  limiited  to  only»  and  can  applied  to  other  minor  atmospheric 
constituents.  The  basic  requirements  for  such  an  application  are  that  spec¬ 
tral  line  features  can  be  isolated,  either  by  a  spectrometer  or  by  a  narrow 
filter,  that  the  corresponding  individual  line  parameters  will  be  known,  and 
that  the  geometry  of  the  experiment  will  define  atmospheric  layers  by  the 
actual  altitudes  or  angles  of  observations.  Line  by  line  -  layer  by  layer 
radiance  analysis  is  then  an  efficient  method  for  deriving  the  vertical  dis¬ 
tribution  of  minor  atmospheric  constituents. 
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VI.  APPENDIX 


It  can  be  shown  that  the  formulation  given  by  Eqs.  (1)  -  (3)  are 
mathematically  equivalent  to  the  more  commonly  used  expression 

R(v,hn)  =Y  B(v,m)  [T(v,  n,  m+1)  -  T(v,n,  m)],  (A-l) 

m=n 

where 


m 

T(v,n,  m)  =  Tf  T(v,j),  (A-2) 

j=n 

and 


T(v,  n,  n+1)  =  1. 

T(v,  n,  m)  is  the  combined  transmittance  through  layer  n  to  m.  Eq.  (A-l) 
is  the  representation  in  terms  of  layers  of  the  integral  form  of  the  radiative 
transfer  equation^^  with  a  zero  boundary  term: 


R(v_hn)=j"  B(v,h)|p^dh,  (A-3) 

n 

where  T(v,h)  is  the  spectral  transmittance  at  height  h. 

Eq.  (3)  is  more  convenient  than  Eq.  (A-l)  for  layer  by  layer  deri¬ 
vations  of  U(n),  while  Eq.  (A-l)  is  more  convenient  when  one  layer  is  as- 
sumed  from  each  altitude  h^  to  the  top  of  the  atmosphere.  However,  in 
both  representations,  improper  evaluation  of  the  transmittance  terms 
might  lead  to  serious  errors  in  U(n).  The  basic  difficulty  is  that  the  ex¬ 
pressions  for  RUhJ  in  Eqs.  (3)  and  (A-l)  are  valid  for  infinite  resolution 
only.  The  infinite  resolution  transmittance  values  for  the  stronger  of  the 
H2^  1*nes  *n  the  ^0  “  30M.m  region  over  atmospheric  paths  of  >  1  km  ap¬ 
proach  the  square  root  region  of  the  curve  of  growth  even  near  30  km  alti¬ 
tude.  This  prohibits  the  use  of  a  number  of  simplifications  that  could  have 
been  made.  For  example,  the  approximation 
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jR{v,hn_1)  t(v,  n)dv~  fR(v,hn_1)dv, 


(A-4) 


might  have  simplified  the  analysis  to 


AR  =  lR(v,  h  )dv  -  fR(v,h  ,  )dv  =  P[1-t(v,  n)]  B(v,  n)dM  (A-5) 

where  AR  represents  the  difference  in  integrated  radiance  observed  at  two 
consecutive  altitudes.  Eq.  (A-5)  reduces  the  problem  to  a  one  layer  calcu¬ 
lation  for  each  layer,  but  the  approximation  in  Eq.  (A-4)  is  not  valid  unless 
T(v^n)  is  in  the  linear  region  of  the  curve  of  growth.  In  order  to  demon¬ 
strate  the  possible  error  associated  with  this  approximation,  trial  runs 
were  made  with  two  layers,  the  first  from  29.  3  km  to  the  top  of  the  atmos¬ 
phere  and  the  second  between  25.  2  km  and  29.  3  km.  The  results  show  that 
U(2)  as  derived  from  Eq.  (A-5)  is  about  4  times  smaller  than  U(2)  as  der¬ 
ived  from  Eq.  (6). 

It  should  be  noted  that  the  square  root  region  implies  a  strong 
pressure  dependence  of  the  emission.  As  a  result,  the  inferred  mixing 
ratio  for  large  enough  layers  will  be  biased  by  the  higher  pressure  levels. 
Trial  runs  with  >10  km  layers  between  12.  5  to  25  km  verify  that  and  yield 
mixing  ratio  values  which  are  close  to  the  minimum  values  obtained  from 
the  10  and  15  layer  calculations. 

Another  difficulty  is  that  the  infinite  resolution  t(v,  n)  cannot  arbi¬ 
trarily  be  replaced  by  some  average  transmittance  values  T(v,n),  as  der¬ 
ived  from  a  band  model  theory  or  from  degraded  spectra.  This  is  due  to 

the  fact  that  R(v,h  )  involves  products  of  the  type  T(v,n,  m)  given  in  Eq. 

n 

(A-2),  and  for  such  products 


! T]"t(v,  i)dv  4=  .f  IT v»  i)dv» 


( A- 6 ) 


unless  all  t(v,  i)  are  in  the  linear  region.  The  error  involved  in  Eq.  (A-6) 

(50-54) 

for  a  two  layer  calculation  has  been  studied,  and  it  has  been  shown 

that  an  overlap  correction  is  needed,  the  neglecting  of  which  can  lead  to 
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large  errors  in  the  inferred  gas  amounts. 

It  follows  from  the  previous  discussion  that  proper  evaluation  of  the 
transmittance  is  critical  for  the  calculations  of  integrated  radiance  for  a 
layered  atmosphere.  In  the  present  analysis,  the  difficulties  have  been 
eliminated  by  performing  exact  line  by  line  transmittance  computations. 
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Fig.  2 
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Table  I 

I 

Balloon  flight  22  February  1971.  Times,  Altitudes,  Pressures 
and  Temperatures  for  selected  records.  The  integrated  radi¬ 
ance  is  for  the  25p,m  line  group. 


Rec 

Time 

MST 

Altitude 

(km) 

Pressure 

(mb) 

Temperature 

(  K) 

Integrated  Radiance 

,  -2  -1! 
(liwcm  sr  ) 

42 

0444 

8.  6 

333.  5 

237.  1 

11.89 

45 

0450 

9.5 

289.0 

229.  5 

5.  03 

47 

0454 

10.  2 

263.0 

224.  0 

4.  36 

49 

0458 

10.  9 

232.4 

217.  6 

3.49 

54 

0508 

12.  5 

180.  0 

216.  5 

2.36 

59 

0517 

13.  9 

141.  0 

216.  5 

1.96 

63 

0525 

15.  2 

116.  6 

216.  5 

1.71 

69 

0537 

17.  0 

88.  5 

216.  5 

1.47 

74 

0546 

18.  3 

71.  9 

216.  5 

1.  30 

80 

0558 

20.  1 

54.  8 

216.  5 

1.  12 

86 

0609 

21.  8 

41.  6 

218.  3 

0.93 

92 

0621 

23.  6 

31.  6 

220.  0 

0.84 

98 

0632 

25.  2 

24.  8 

221.  6 

0.  66 

104 

0644 

26.  9 

18.  9 

223.  3 

0.  57 

127 

0728 

29.  3 

13.4 

225.  6 

0.  52 

Table  II 


Balloon  flight  29  June  1971.  Times, 
Temperatures  for  selected  records, 
is  for  the  25ym  line  group. 


Rec 

Time 

Altitude 

Pressure 

MST 

(km) 

(mb) 

28 

0228 

9.3 

316.  0 

29 

0230 

9.8 

291.7 

30 

0232 

10.4 

269.4 

31 

0234 

11.  0 

247.  6 

37 

0246 

14.4 

143.  5 

40 

0252 

15.9 

113.  2 

43 

0257 

17.  1 

92,  0 

48 

0307 

18.  9 

68.  3 

51 

0313 

20.  1 

56.  3 

54 

0319 

21.  3 

46.  7 

58 

0326 

22.  8 

37.  1 

61 

0332 

24.  0 

30.  8 

65 

0340 

25.  6 

24.  2 

69 

0348 

27.  1 

19.  2 

74 

0357 

28.  8 

14.  8 

31 


Altitudes,  Pressures  and 
The  integrated  radiance 


Temperature 


Integrated  Radiance 

.  -2  -1 
(pwcm  sr  ) 


241.  2 
237.0 
232.  7 
228.  3 

205.8 

204.  7 

205.  6 
208.  9 
211.  3 
214.  9 
218.4 
220.  2 
222.  6 

224.9 
227.  5 


11.  14 
9.  12 
7.  19 
5.  60 
1.46 
1.22 
1.15 
1.  06 
1.  00 
0.  91 
0.79 
0.73 
0.  59 
0.  55 
0.48 


